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Tracing of Skew Rays 

By Donald P. Feder 

This paper presents a set of convenient recurrence formulas for the tracing of skew rays 
by means of a desk calculator. Because of the increased speed and field coverage of modern 
lenses, it is frequently desirable to trace such rays in making a complete study of an optical 
system. 

The formulas apply to systems of spherical refracting surfaces. The system need not 
be centered. 

Particularly convenient check formulas are available, which need be applied only at 
the final surface. 



A complete study of an optical system appears 
to demand some knowledge of the behavior of the 
general or skew rays. A computational process 
for obtaining this information by tracing such rays, 
which has proved convenient in practice, will now 
be described. 

It is contemplated that these formulas will be 
used with a desk calculating machine. When such 
a calculator is used, twelve numbers per surface 
must be written down for each ray. 

An advantage of the method is the existence of 
a simple check that need be applied only at the 
final surface. 

Brief derivation. We consider a system of 
spherical refracting surfaces. The system need 
not be centered. 

We denote the incident and refracted rays by 
the unit vectors Q and Q', and the radius vector 
to the center of curvature from the point of 
incidence by R. Bold-faced type will be used to 
denote vectors. 

We define a new vector P 



P=N QXR. 

The law of refraction becomes 
P=N' Q'XR, 



(1) 



(2) 



where N and N' are the refractive indices of the 
surrounding media. 

We can now choose our independent variables 
from the set of vectors (P, Q, R). Any two of 



these determine the direction and location of the 
ray and may be used as independent variables. 
An obvious choice is the pair Q and R, but P and Q 
yield simpler equations, particularly a simpler 
transfer. 

We eliminate R from eq 1 and 2 in order to obtain 
Q f as a function of P and Q. Multiplying eq 1 vec- 
torially by Q, one can obtain R, which is then 
substituted in eq 2. After solving the resulting 
equation for Q' 7 we obtain 



where 



Q'=aQ+bPXQ, 

P 2 +NN'(Q'.R)(Q.R) 

NN'R 2 

N'(Q'.R)-NjQ-R) 

NN'R 2 



(3) 
(4) 
(5) 



(Q'R) and (Q'-iJ) may be evaluated by squaring 
eq 1 and 2 respectively. We find 



N 2 (Q-R) 2 =N 2 R 2 -P 2 



N'\Q'-R) 2 =N' 2 R 2 - 



(6) 

(7) 



Given P and Q, these equations enable us to 
find Q' = Qi. It remains only to find Pi for the 
following surface in terms of P and Q at the 
preceding one. 

In figure 1 and ft are the respective centers, 
and A and A x are the points of incidence of the 
ray. In general the quadrilateral is not plane. 
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P\x — Px 



Figure 1. 

R, Q\ 9 and R\ denote vectors. 



Denoting the vector from to X by T, one 
has immediately 



Ri+AAiQi = R+T 



(8) 



Multiplying vectorially by N'Q'=NiQ h the term 
in AA\ vanishes, and we have 



Pi=P+QiXt 



(9) 



where t=N 1 T 



Summary in scalar form. Since the vast ma- 
jority of optical systems are centered, let lis sup- 
pose that 00i lies on the optical axis. We may 
take the x axis along this line so that t y =t z =0. 

Suppose that Q=(X, Y, Z) and P = (P X , P V) P z ) 
are given in scalar form. Defining the auxiliary 
constants r=N\R\ and r'=N'\R\, and writing 



P 2 =P 2 +P 2 + p2 

g *=r 2 -P 2 
g*=r' 2 -F- 

„_g'g+P 2 



b= 



rr 
9' -9 



rr 



X^aX+blPyZ-PtY) 
Y l =aY+b[P,X-P x Z] 

Z 1 =aZ+b\P x Y-P v X] 
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(10) 

(ID 
(12) 

(13) 
(14) 

(15) 



Ply Py-\-t x Zl 

Plz — Pz t x Yi 



(16) 



The signs of g and g' are identical and are posi- 
tive for surfaces convex toward the entering rays. 

Equations 10 through 16 inclusive, form the 
complete set for systems of spherical surfaces. 

One notices that P x is an invariant of the ray. 
P x =0 implies that the ray lies in a meridian plane. 

Plane surfaces. At a plane refracting surface 
normal to the optical axis, the refraction equa- 
tions are greatly simplified. 

From eq 1 and 2 by taking R = i one obtains 



also 



N'Y'=NY 
N'Z'=NZ 

2 2 

X' =1-F' -Z' 



(17) 
(18) 



and X f is always positive. 

Since P does not define the location of the ray 
at a plane, we introduce a vector F, lying in the 
plane and extending from the axis to the inter- 
section of plane and ray. 




Figure 2. 
R, Qi, F, and T denote vectors. 

From figure 2 we obtain AAiQ x = R+T+F and 
multiplying by NiQ x as before, we have 

P^NxFXQi+tXQi (19) 

where t=N\T. In scalar form, if we take t y =t 2 =0 

P^NAFyZ-FzY,] n 
P y =N 1 F z X 1 -t x Z 1 V (20) 

P 2 =-N 1 F y X 1 +t x Y^ 
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Evidently eq 20 may be used to obtain F in terms 
of P or conversely. Frequently they are used as 
initial and final equations as well as at refracting 
surfaces. 

Checking equations. At each surface 






(21) 



from eq 1 or 2. It is important to note that these 

checks need be applied only at the final surface. 

For we have from eq 3 by scalar multiplication 

PQ'=aPQ 

and similarly from eq 9 

/V0i = P-Qi=aP-0. (22) 

Furthermore, since a can be shown to be always 
nonzero, it follows that if PQ is nonzero at any 
point, it remains nonzero. This implies that an 



error that causes PQ to differ from zero at any 
surface will appear at the final surface. 

Similarly, if any error is made so that Q 1 — 1^0, 
the error will appear at the final surface. 

In addition one must have at each surface 



a 2 +b 2 P 2 -l = 



(23) 



These checks should guarantee the accuracy of 
the final data, provided no errors are made in the 
constants. 

Using this method, the accuracy of the results 
is not limited by a particular set of tables. All 
that is needed is two square roots, and this may 
be obtained, with any desired accuracy, by the 
usual iteration formulas. It is very easy to ob- 
tain results to 10 significant figures by employing 
at most a small table. 

Washington, March 30, 1950. 
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